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In a Stark decelerator, polar molecules are slowed down and focussed by an inhomogeneous electric 
field which switches between two configurations. For the decelerator to work, it is essential that 
the molecules follow the changing electric field adiabatically. When the decelerator switches from 
one configuration to the other, the electric field changes in magnitude and direction, and this can 
cause molecules to change state. In places where the field is weak, the rotation of the electric 
field vector during the switch may be too rapid for the molecules to maintain their orientation 
relative to the field. Molecules that are at these places when the field switches may be lost from the 
decelerator as they are transferred into states that are not focussed. We calculate the probability of 
nonadiabatic transitions as a function of position in the periodic decelerator structure and find that 
for the decelerated group of molecules the loss is typically small, while for the un-decelerated group 
of molecules the loss can be very high. This loss can be eliminated using a bias field to ensure that 
the electric field magnitude is always large enough. We demonstrate our findings by comparing the 
results of experiments and simulations for the Stark deceleration of LiH and CaF molecules. We 
present a simple method for calculating the transition probabilities which can easily be applied to 
other molecules of interest. 

PACS numbers: 37.10.Mn, 37.20. +j, 31.50.Gh 



I. INTRODUCTION 

Beams of fast-moving polar molecules can be brought 
to rest using a Stark decelerator [1^. The molecules ex- 
perience a force in the inhomogeneous electric field of 
the decelerator. By switching the field on and off in a 
carefully tailored sequence this force always opposes their 
forward motion. The Stark decelerator was first demon- 
strated for metastable CO [2] and several other molecules 
have since been decelerated this way. The decelerator 
works best for molecules in low-field seeking states since 
they are easy to focus through the device. For the de- 
celerator to work it is crucial that the molecules remain 
in the same quantum state throughout the deceleration 
process. The timing sequence that needs to be applied to 
the decelerator electrodes depends on the Stark shift and, 
since different states have different shifts, a molecule that 
changes state will no longer be decelerated. Moreover, if 
the state change is to one of the high-field seeking states, 
the molecule will be attracted to the electrodes where the 
field is strongest, and it is unlikely to be transmitted. 

The unwanted change of state can occur because the 
molecules experience a time- varying electric field as they 
move through the decelerator and as the decelerator is 
switched from one field configuration to another. Tran- 
sitions between states will only be significant if the 
time- varying electric field has frequency components that 
approach the frequency splitting between those states. 
Some simple order-of-magnitude refiections provide con- 
siderable comfort in this respect - in a typical deceleration 
experiment the Stark shifts are tens of GHz, whereas the 
molecule takes about 10 /iS to travel from one stage to the 
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next and the rise and fall times (t) of the switching need 
only be 10 times smaller than that for the deceleration 
to be efficient. Nevertheless, we must consider carefully 
those cases where energy levels happen to be very close, 
as occurs at an avoided crossing, or when the electric field 
becomes so small that states differing only in their ori- 
entation with respect to the field are nearly degenerate. 
It is important to distinguish these two cases. In the de- 
celerator, molecules need to traverse an avoided crossing 
diabatically so that the gradient of the Stark shift is the 
same on both sides of the crossing. That may be difficult 
to arrange if the energy gap at the avoided crossing is too 
large. Conversely, to preserve its orientation in low field, 
the molecule must follow the changing field adiabatically, 
which can be arranged by ensuring that the electric field 
in the decelerator is never too small. 

When the decelerator switches, both the magnitude 
and the orientation of the electric field changes. It is the 
rotation of the electric field that is responsible for the 
orientation-changing transitions. On the axis of the de- 
celerator the electric field vector rotates through 90° each 
time the decelerator switches. The rotation rate of the 
field vector depends on the position within the decelera- 
tor. At positions where the field magnitude is small prior 
to switching and large afterwards (or vice versa) , the ma- 
jority of the 90° rotation occurs in a small fraction of r, 
meaning that the rotation rate reaches values very much 
larger than r~^. In fact, as we shall see, the rotation rate 
is at its largest at a point in space and time where the 
energy splittings are at their smallest and the molecules 
most vulnerable to nonadiabatic transitions. It is these 
rotation-driven transitions that are the main subject of 
this paper. 

Losses due to nonadiabatic transitions and due to 
avoided crossings have been considered previously in the 
literature on Stark deceleration. For the deceleration of 
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OH and NH in low-field seeking states, there is good ex- 
perimental evidence that nonadiabatic transitions do not 
occur [21 13]. For the specific case of OH, this is supported 
by theoretical work However, for LiH in the first ro- 
tationally excited state, numerical simulations could only 
be reconciled with experimental results by including in 
the simulations the probability of nonadiabatic transi- 
tions from the low-field to the high-field seeking compo- 
nent, driven by the rotation of the electric field during 
the switch [6] . In experiments on the deceleration of OH 
molecules in high-field seeking states, there was strong 
loss due to a crossing of hyperfine states which was elim- 
inated by applying a small bias voltage so that the de- 
celerator field was always above the field where the two 
states cross ^7^. In experiments on SO2, potential loss 
due to an avoided crossing of the decelerated state was 
considered, but occurred at too high a field to be im- 
portant [8j. Loss at an avoided crossing has been dis- 
cussed in the context of decelerating asymmetric rotors 
[S] . Loss of molecules due to nonadiabtic transitions was 
also found to be significant in the recently demonstrated 
chip decelerator |10j . and is important in traps for polar 
molecules, as observed in [TT] and discussed theoretically 
in 0113 El]. 

We will first outline the theory of nonadiabatic transi- 
tions in a Stark decelerator. This leads to a simple adia- 
baticity criterion and to a simple method for calculating 
transition probabilities. We then provide some example 
calculations, showing that the transition probability is a 
strong function of position in the decelerator, and that 
molecules in the non-decelerated group are far more likely 
to be lost by this mechanism than those in the deceler- 
ated group. Finally, we compare the theoretical results 
with the results of experiments on the Stark deceleration 
of LiH and CaF molecules. 



II. THEORY 

Figure [T]ja) illustrates the layout of a typical deceler- 
ator P], showing three deceleration stages. Each stage 
is formed from a pair of parallel rods with the molecu- 
lar beam passing between them. The axes of the rods 
are parallel to the z-axis in the odd-numbered stages 
and parallel to the x-axis in the even-numbered stages. 
The molecular beam travels along the y-axis and, to ex- 
press the periodicity of the decelerator, the longitudinal 
coordinate is written (in radians) as the reduced posi- 
tion 6 = TTjj/L, L being the distance between successive 
stages. The decelerator is switched between two states, 
state 1 which has the odd stages grounded and the even 
stages charged (as shown in Fig.jlja)), and state 2 which 
has the even stages grounded and the odd stages charged. 
The profile of electric field strength in these two states is 
shown in Fig.[l]jb) for a decelerator that has i? = 1.5 mm, 
ro = 1 mm, L = 6 mm and V ^ 15 kV. On the y-axis, the 
electric field is directed along z in state 1 and along x 
in state 2. When the decelerator switches between these 



two states the field on the beamline rotates by 90° and 
changes in magnitude from the value given by the solid 
line to that given by the dashed line in Fig.[l|b). 

The switching sequence applied to the decelerator is 
synchronized to the motion of a molecule having a par- 
ticular initial position and velocity. By construction, this 
synchronous molecule has the same value of 9 (modulo 
tt) every time the decelerator switches, and that value is 
called the synchronous phase angle, (/>o. In the coordi- 
nates of Fig.[l] (f)o = corresponds to zero deceleration, 
and 00 — 90° corresponds to maximum deceleration. The 
decelerator exhibits phase stability jT3], meaning that it 
traps a group of molecules centred in phase space on the 
synchronous molecule. These phase-stable molecules re- 
main close to the synchronous molecule throughout the 
deceleration process. The area they occupy in phase 
space is the longitudinal acceptance of the decelerator 
and depends on the degree of deceleration applied. Fig- 
ure [TJc) shows an example of this acceptance for several 
values of (pQ. We will return to this figure to explain how 
the phase-stable molecules are less prone to nonadiabatic 
transitions than the non-phase-stable molecules. 

We consider a rigidly rotating linear molecule of dipole 
moment /x and rotational constant B in an electric field E 
which is Eiz before the switch and i?2X after the switch. 
During the switch, the first field decays with a time con- 
stant Ti while the second field turns on with time con- 
stant T2. The time-dependent Hamiltonian, H{t), that 
describes this system is given by 

H{t) = hBN^-fiE,e-*/-'C^a'^-^{l~e-'/^-){c['l^ci'l), 

v2 

where N is the rotational angular momentum operator, 
and we have written the angular factors in the Hamilto- 
nian in terms of the components of the rank 1 spherical 

harmonic operator, Cq^^ = \J^Yq'^\ 

In the absence of the electric field, the eigenstates 
are the angular momentum eigenstates \N^M), satis- 
fying N^\N,M) = N{N + l)\N,M) and N^\N,M) ^ 
M\N^M). The electric field mixes states of the same 
M but different N so that N is no longer a constant of 
the motion. We use the notation jA/", M) e to denote the 
mixed- eigenstates in the static field Ei, omitting the 
subscript E when it is not needed. To identify a specific 
state we assign a value to N which is equal to TV in the 
limit that i? — >■ adiabatically. These eigenstates, and 
the corresponding eigenvalues, ea/'.m, are calculated by 
diagonalizing -ff(O) for many different values of Ei us- 
ing a suitably large basis set formed from the field-free 
eigenstates. Figure [2] shows the energies of the lowest- 
lying states as a function of the electric field. Both the 
energy and the electric field are written in dimensionless 
units by dividing the former by the rotational constant 
hB, and the latter by the ratio of rotational constant to 
dipole moment hB/fi. The Stark decelerator works best 
for the states with the largest positive Stark shift since 
these states are naturally focussed and are most strongly 



3 




Reduced position, 6 (degrees) 





Reduced position, 6 (degrees) 

FIG. 1: (Color online) (a) Structure of the Stark decelerator, 
showing three consecutive deceleration stages centred on the 
n'th stage. Each stage is formed from a pair of rods whose 
axes are parallel to the z-axis for the odd-numbered stages, 
and to the a;-axis for the even- numbered stages (n is even). 
For clarity, the two rods forming the n'th stage are shown 
slightly displaced, though in reality one should be exactly 
behind the other, (b) Electric field profile versus reduced po- 
sition for the two switch configurations when R — 1.5 mm, 
ro = 1mm, L = 6 mm, V = 15 kV. (c) Longitudinal phase- 
space acceptance for synchronous phase angles of 0° (outer- 
most line), through to 80° (innermost line), in steps of 10°. 
The plot is for CaF molecules in the |4, 0) state and for the 
electric field profile shown in (b). 



decelerated. Note that these low-field seeking states all 
turn over to become high-field seekers when the field is 
large, and that the turning point shifts to higher electric 
field for larger values of Af and B. 

At time t — the molecule is in the particular eigen- 
state \J\fp, Mp) We can calculate the time evolu- 
tion of the state vector \a) in the z-fixed basis \NM), 




Electric field /(hB/n) 



FIG. 2: Stark shift of the low-lying states of a rigid rotor 
molecule. States are labelled by their values of A/" and M. 



by solving the Schrodinger equation for the coefficients 

CNM = {NM\a), 

ihjCNM= Y.^^^"^\H{iWM')cN'M', (2) 

N'M' 

with this initial condition and with H{t) given by Eq. ([l]). 
Having found the time evolution of \a), we can then cal- 
culate the probability of finding the molecule in any one 
of the states \N , M) ^(t) as a function of t. 

Although this calculation is straightforward, the in- 
tegration of Eq. ([2]) is slow because the basis set needs 
to include a large number of N states and because the 
timescale of the switch is vastly larger than \/B, which 
sets the scale for the time step in the integration. For- 
tunately, a few simple transformations provide us with 
an equation that is much faster to solve and which also 
makes it obvious that it is the rotation of the electric field 
vector that drives the transitions. 

If the coordinate system were rotated so that the field 
at time t pointed along the local z'-axis, the instanta- 
neous Hamiltonian would be H' = hBN'^ — fj,E{t)CQ^\ 
E{t) being the electric field magnitude at time t. Hamil- 
tonian ([1]) is related to this one by a rotation about the 
y-axis through an angle /3(t), H = D~^H'D. Here D 
is the relevant rotation operator, D = exp{—iNy/3/ti). 
The evolution of the state vector \a) is governed by 
the Schrodinger equation: ih^\a) = H{t)\a). Apply- 
ing the rotation operator to both sides of this equation, 
and introducing the state vector in the rotating frame, 
I a') — D\a), we obtain 

ihD^{D-^\a')) = DH\a) = bb~^H'b\a) = H'\a). 
* (3) 

Using 

dt h y dt ' ^ ' 

and the fact that commutes with Ny, we obtain the 
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time evolution of the state in the rotating frame, 



d 



H' +Lu{t)Nyj\a') ^ih-\a'), (5) 



_ d/3 



where a;(t) — . 

Next, following the algebra set out in App.jXJ we ex- 
press this equation in a basis formed from the instanta- 
neous eigenvectors of H' , obtaining 



Uj{t)Ny 



dH' 



dt 



ih 



dt 



(6) 

where hhJmm' — £»« ~ Cm', a„i — (TO|a')e"™*'''', m stands 
for the pair of quantum numbers A/", M, and m' for 
AA', M'. 

In Eq. (|6| the dynamics are expressed in a particularly 
transparent way. The first operator in the matrix element 
is due to the rotation of the electric field vector and cou- 
ples states of different M, while the second arises from 
the change in the electric field magnitude and can only 
couple states of the same M belonging to different man- 
ifolds Af. Provided the matrix elements in the sum vary 
little on the timescale of l/uj^m', the terms in the sum 
oscillate around zero with an angular frequency Umm'- 
Those amplitudes am that are initially zero will always 
remain small unless they grow significantly in a time of 
order l/wmm'- Thus, if the matrix elements are all much 
smaller than the energy splittings hojmm', tbe state am- 
plitudes hardly change meaning that the molecule adi- 
abatically follows the changing field. This is the situa- 
tion that we want in the Stark decelerator, so let us now 
estimate the transition probability arising from the off- 
diagonal terms. Assuming that these terms are at least 
small, we can treat them using time-dependent first-order 
perturbation theory. For a perturbation V, the probabil- 
ity Pn of being in state |n), having started out in the 
initial state is at most 4| (n|F|i) ^/(fia;™)^ (see [15] 
for example). 

Taking first the last operator on the right hand side 
of Eq. (jei, we have {n\V\i) = ifi'§{n\C^^'>\i)/uJm. The 
maximum value of this last matrix element is 1, while 
that of ^ is -Emax/T, i^max being the maximum electric 
field in the decelerator. As shown in Fig.[2j the weak- 
field seeking states have turning points, and this sets an 
upper limit to i?max- Within manifold A/", the state with 
the largest positive Stark shift is \Af, M = 0) , whose 
turning point is at a field of approximately ATrhBJ\f^ / ^. 
For those states that are coupled by the perturbation, the 
smallest possible value of Wni occurs at low electric field 
and is AirBAf. Putting all this together, we conclude that 
the probability of a transition induced by the changing 
electric field magnitude is smaller than l/(27ri?T)^. Most 
molecules of interest have B > 10 GHz, while the rise 
and fall time of the switch is typically t > 100 ns, and 
so the maximum transition probability, about 10^*^, is 
completely negligible. 



Now we turn to the rotation of the field, i.e. to the 
matrix elements of ujNy in Eq. Here, the situation 
is quite different because this term couples states of dif- 
ferent M within the same manifold, and these states ap- 
proach each other very closely when the field is small. 
The matrix element that we need is only non-zero when 
M' = M ± 1, in which case it is 

(A/-, M±l\Ny\N', M) = ^ ^{NtM){N±M + 1) 

N 

{M,M ±l\N,M ±1){N,M\N,M). (7) 

At low electric field, where transitions are most likely to 
occur, the mixing of different N states is small and so 
only one term {N = M) makes any significant contribu- 
tion to the sum. The probability of finding the molecule 
in the state M' = M ± 1 is therefore 



P±i(t) < -^{UtM){U±M + 1). 



(8) 



For low-lying rotational states, which are the ones of most 
interest to us, the transition probability will always be 
very small if a; ^ Uni - this is the criterion for the rota- 
tion to be adiabatic. 



A. Examples 

We consider two examples, LiH prepared in the ini- 
tial state \i) = |1,0), and CaF prepared in the initial 
state \i) — |4,0). They move through the decelerator 
illustrated in Fig.[IJa), with R = 1.5mm, ro = 1mm, 
L = 6mm and y = 15 kV. We take \n) = |1, 1) for LiH 
and \n) = |4, 1) for CaF, and we take the rise and fall 
time of the switch to be r = 500 ns. Figure [3] shows 
how the values of a; = d(i/dt (solid lines), (jj}^f^ (dotted 
lines) and w^f^ (dashed lines) change as a function of 
time for several different values of the reduced position 



When e = -90°, cj exceeds 



,CaF 



for the first 15 ns 



after the switch is thrown. At this position in the decel- 
erator the field changes from a very small field along z 
to a very large field along x (see Fig.[IJa,b)), and so the 
field vector completes almost all of the 90° rotation in 
the first few ns, i.e. in a small fraction of r. As a result, 
the rotation rate is very large over this interval of time, 
Lu ^ 10* rads~^. At the same time, the total field is still 
very small, and so the angular frequency interval between 
the two M states is also very small, w,^f^ ~ lO^rads"^. 
We see that the rotation is far from adiabatic. A CaF 
molecule that happens to be at this position when the 
field switches is very likely to end up in a different state. 
For LiH the situation is not as severe because, at early 
times, iJ^f^ exceeds w by a factor of 3. Nevertheless, 
since the two frequencies are still of the same order of 
magnitude, we can expect there to be a significant tran- 
sition probability. As we move away from this position 
in the decelerator, uj decreases as the difference in field 
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FIG. 3: (Color online) The rotation rate, dfi/dt (solid), the 
angular frequency separation of LiH |1,0) and |1,1) states 
(dotted), and the angular frequency separation of CaF |4, 0) 
and 1 4, 1) states (dashed), as a function of time since throwing 
the switch and at various reduced positions (8) within the 
decelerator. The switch rise and fall time was set to r = 
500 ns, and the electric field profile used in the calculation 
is the one shown in Fig.[l](b). Note the logarithmic scale on 
both axes. 



magnitude between the two switch states decreases, and 
LOni increases because the minimum field value increases. 
At 9 = —54° for example, the maximum value of lo is 
about 3 times smaller than the minimum value of w^f^, 
while at 9 = —18° the former is more than two orders 
of magnitude smaller than the latter. At all positions, 
the Stark splitting for LiH is far greater than for CaF 
meaning that LiH is not as susceptible to nonadiabatic 
transitions. The situation at 6' = 90° again seems dan- 
gerous - the field starts out very large along z and ends 
up very small along x. This means that most of the rota- 
tion occurs at the very end of the switch time, at which 
time the frequency separation of the two states is very 
small. However, because the gradient of the exponen- 
tial rise and fall is far smaller at late times than at early 
times, the rotation is not as fast and the maximum value 
of cj/cjCaF is about 1/20. 

To calculate the exact time-evolution of the state of 
the molecule, we set Af' — Af in Eq. (6]), having already 
established that couplings between different Af manifolds 
are negligible. Then we solve numerically the resulting 
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FIG. 4: (Color online) Evolution of the populations in the 
sublevels of the A/" = 1 manifold of LiH as a function of time 
since throwing the switch and at the position = —90°. The 
curves are labelled by their value of \M\. The initial state is 
M — 0, the switch rise and fall time is r = 500 ns, and the 
electric field profile is the one shown in Fig.[T][b). 



set of differential equations, 



AT 



(TV, M'\uj{t)Ny\Af, M)aMe"''"(-^-"'<'^'^'''* = ih 



daM' 
dt 

(9) 

Figure |4] shows what happens to LiH molecules ini- 
tially in the |1,0) state, when they are a.t 9 — —90° at 
the moment of switching. Here, the rise and fall times 
of the switch are 500 ns, and the electric field profile is 
as shown in Fig.ljb). As anticipated from the consid- 
erations above, the evolution is not fully adiabatic. The 
rapid rotation of the field vector that occurs during the 
first 10 ns results in some population transfer into the 
A/ = ±1 states. After a few tens of ns the evolution 
becomes adiabatic and there is no further change in the 
state amplitudes. Moving away from this 9 — —90° posi- 
tion we find that the transition probability drops rapidly 
to zero. We note that in all cases the amplitudes of the 
±M states are identical because the molecule starts in 
M = and, in our coordinate system, the z-axis lies 
in the plane of rotation. This symmetry is also seen by 
setting Af = in Eqs. ([7| and ([8|. 

Figure [5] shows what happens to CaF molecules ini- 
tially in the |4, 0) state, when they are at 6* = —90°, 
-70° and -50°. When 9 = -90°, the initial popula- 
tion in the |4, 0) state is rapidly transferred to the other 
M states within Af = 4. First, the \M\ = 1 state is 
populated, and from here the \M\ — 2,3 and 4 states 
all become populated as time progresses. After about 
50 ns the populations stop changing because the evolu- 
tion changes over from being strongly nonadiabatic to 
being adiabatic, as anticipated from our consideration of 
Fig.jsja). At 6' = —70° the evolution is slower and less 
of the population is driven to the higher \M\ states. In 
particular, the \M\ — 3 states, which received the largest 
share of the population aX 9 — —90°, now receive very 
little, and the \M\ = 4 states are never populated at all. 
These trends continue as 9 increases, the evolution get- 
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FIG. 5: (Color online) Evolution of the populations in the 
various sublevels of the A/" = 4 manifold of CaF as a function 
of time since throwing the switch, and at three different re- 
duced positions in the decelerator. The curves are labelled by 
their value of \M\. The initial state is M = 0, the switch rise 
and fall time is r = 500 ns, and the electric field profile is the 
one shown in Fig.[T]^b). 



ting progressively slower and the population being con- 
fined more and more to the low \M\ states. At 9 = —50° 
all the population remains in the M = and \M\ = 1 
states. Again, this is consistent with our expectations 
from Fig.jsjb) and Eq. Q. It is also interesting to study 
the situation at 6* = 90° , where it is necessary to find the 
solution up to i ^ 3 /iS so as to capture any dynamics that 
may occur when uj and ujni pass through their point of 
closest approach (see Fig.jsjjf)). We find the population 
transfer to be negligible here. 

Figure |6] summarizes the results we have obtained. It 
shows the transition probability for LiH |1,0) and CaF 
1 4, 0) molecules as a function of the position in the decel- 
erator at the moment of switching. Transitions occur in a 
zone centred on 9 = —90°. For CaF the maximum tran- 
sition probability is 84% and the zone has a full width 




%0 



-45 45 

Reduced position, 9 (degrees) 



FIG. 6: Transition probability out of the M — state as a 
function of the reduced position in the decelerator. The po- 
sition dependence is symmetric about 9 — —90°. The switch 
rise and fall time is r = 500 ns, and the electric field profile is 
the one shown in Fig.jljb). Solid line: LiH in A/" = 1. Dashed 
line: CaF in A/" = 4. 



at half maximum of about 60°. For LiH, the maximum 
probability is 11% and the zone has a full width at half 
maximum of about 20°. Obviously, faster switching or 
lower applied voltages will tend to increase the transi- 
tion probabilities and extend the width of the zone. 



B. Loss 

In discussing the degree of molecule loss to be ex- 
pected from these unwanted transitions, we need to dis- 
tinguish molecules that are phase-stable from those that 
are not. The phase-stable molecules lie within the longi- 
tudinal phase-space acceptance of the decelerator; they 
are the ones being decelerated. The phase angle of these 
molecules, called 4> and meaning the value of 9 at the mo- 
ment of switching, is constrained by the action of the de- 
celerator to lie within the boundaries of the phase-space 
acceptance. These boundaries are shown in Fig.[lJ^c) for 
various values of the synchronous phase angle, (jjQ. This 
figure shows that if the decelerator is operated at large 
phase angles, say 00 > 45°, corresponding to strong de- 
celeration, the phase-stable molecules cannot reach the 
region centred on 9 = —90° where transitions are likely 
to occur. The group of molecules being decelerated are 
largely protected from transitions driven by the switch- 
ing field. If we plot a molecule's phase-space position 
every time the decelerator switches we find that the 
phase-stable molecules move on closed trajectories cen- 
tred on the synchronous molecule in the phase-space plot 
of Fig.jljc). Thus, even when the decelerator is operated 
at low phase angle, there will be a sizeable fraction of 
phase-stable molecules that never enter the regions where 
transitions are likely. By contrast, the molecules that are 
not phase-stable tend to explore all possible values of (j) 
as they make their way through the decelerator. Since 
the region where transitions occur is sizeable, and since 
in a typical deceleration experiment the field switches 
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FIG. 7: The Stark shift of the energy levels of CaF that cor- 
relate to A*' = 4 at low field. The states are identified by their 
values of J and F at low field, and by the value of |Mj?|. At 
higher field the states separate into 5 groups that can be la- 
belled by the projection quantum number |A/]v|. Those that 
join the Mjv = group are shown by bold black lines. The 
inset shows the same set of states up to a field of 200kV/cm. 

about 100 times, we can expect a large fraction of those 
molecules to be driven into another state and be lost from 
the decelerator. So we reach the interesting conclusion 
that the loss will be far greater for the non-decelerated 
group of molecules than for the decelerated group. 

C. Hyperflne structure 

So far in our discussion, we have neglected hyperfine 
structure. This is acceptable when the hyperfine struc- 
ture is small compared to the Stark shift at all relevant 
fields, as is the case for ''LiH in its electronic ground 
state. Here, the hyperfine splittings are of order 100 kHz, 
and arise from the interaction of the electrons with the 
electric quadrupole moment of the Li nucleus, and from 
the interactions between the nuclear magnetic dipole mo- 
ments and the magnetic dipole moment associated with 
the rotation. These splittings are small enough to be 
neglected in our calculations of nonadiabatic transitions. 
However, in CaF and many similar molecules, the hy- 
perfine structure is comparable to the Stark shift at the 
fields strengths where we are troubled by nonadiabatic 
transitions. This complicates the matter quite consider- 
ably 

In CaF the main hyperfine interactions are a spin- 
rotation interaction, which splits a rotational state N 
into two states labelled by J = N±l/2, and the interac- 
tion between the electron spin and the spin of the fluorine 
nucleus, which further splits each state J into two, la- 
belled by the total angular momentum quantum number 
F — J±l/2. Each state F has 2F + 1 magnetic sublevels 
labelled by Mp, which is a good quantum number at all 
fields. In the absence of a magnetic field, states differing 



only by the sign of Mp are degenerate. Figure [7] shows all 
the hyperfine states correlating to = 4 at low field, as a 
function of electric field. The main figure shows the Stark 
shifts at low fields, while the inset shows the same states 
up to a field of 200kV/cm. The energies were calculated 
by adding to H{0) (Eq. (jl])) the usual hyperfine terms for 
a^E state, iJhyp = 7S-N-h6I-S-hc(i-k)(S-k)-|-CI-N, 
where S and I are the electron and nuclear spin angular 
momentum operators, k is a unit vector along the inter- 
nuclear axis, and 7, b, c and C are constants given for 
CaF in 16]. 

At very low field the energy levels group into 4 well 
separated components labelled by J and F because the 
Stark shift is small compared to the hyperfine structure. 
At high fields the energy levels regroup into 5 well sepa- 
rated components, labelled by the value of |Mjv|, because 
the Stark shift is far greater than the hyperfine splitting. 
The tangle of levels at intermediate fields complicates our 
calculation of the transition probabilities. Instead of pre- 
senting this calculation, we make some general observa- 
tions. First let us concentrate on a single (J, F) manifold. 
It splits up into Mp states in much the same way as do 
the Mn states when we ignore the hyperfine structure 
(i.e. as in the inset). The level ordering is the same, and 
at any given field (over the range shown in the plot) the 
sizes of the splittings between the M states are similar 
to those obtained in the absence of hyperfine structure. 
Since the matrix elements for the transitions between the 
M states are also similar in the two cases we can conclude 
that the transition probabilities between Mp states be- 
longing to the same {J,F) are not very different from 
the Af-changing probabilites that we have already calcu- 
lated. Continuing to focus on just one {J,F) manifold, 
we note that a change of Mp can change the character 
of the state from a low-field seeking to a high-field seek- 
ing one (compare the main figure with the inset). Again, 
this is no different from the situation with the hyperfine 
structure neglected. These considerations lead us to rea- 
son that loss due to nonadiabatic transitions is at least as 
severe when the hyperfine structure is included as when 
it is neglected. In addition to state changes within a 
{ J,F) manifold, transitions might occur between states 
from different manifolds, if these approach or cross one 
another. However, there is only a coupling if the values 
of Mp differ by one. Focussing on the states with the 
largest positive Stark shift (bold black lines, all corre- 
lating to Mn = at high field), which are the ones of 
most interest to us, we see that such crossings or near ap- 
proaches never occur. The same is true of the states that 
correlate to |Afjv| = 1. It is not true of the other states, 
but these are lost in the decelerator anyway. Putting 
these observations together, we conclude that our calcu- 
lations of transition probabilities remain good estimates 
when the hyperfine structure is included. In most cases 
such rough estimates are sufficient since our task is usu- 
ally to understand how severe the loss is likely to be and 
what can be done to avoid that loss. 
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III. EXPERIMENT 

In this section, we discuss experiments on the Stark 
deceleration of ^LiH in the |1, 0) state and of CaF in the 
|4,0) state, in order to demonstrate our main findings 
regarding nonadiabatic transitions. 

A. Setup 

The experimental setup for producing, decelerating 
and detecting LiH has been described before [H |T7] . We 
have not reported the deceleration of CaF previously, but 
the decelerator we use has the same structure as in [6] 
and the methods we use to produce and detect CaF are 
described in [T5|. Both molecular species are produced 
by laser ablation of a target into a supersonically ex- 
panding gas jet. For the LiH experiments, the target is 
Li and the carrier gas is a mixture of H2 (« 2%) and 
Kr. For the CaF experiments the target is Ca and the 
carrier gas is a mixture of SFg (sa 2%) and Xe. The 
temperature and pressure in the solenoid valve are ap- 
proximately 300 K and 3-4 bar respectively. The valve 
is fired at 10 Hz and emits gas pulses with a duration of 
approximately 200 fis into the source chamber which is 
maintained below 10~^ mbar. The fundamental output 
of a Q-switched Nd:YAG laser, with a pulse duration of 
10 ns and an energy of 25 mJ, is used to ablate the tar- 
get outside the nozzle of the valve. The resulting beam 
passes through a 2 mm diameter skimmer and into a sec- 
ond chamber which houses the decelerator and the detec- 
tor. In the LiH source, most of the molecules are formed 
in the ground state and so we excite the molecules to the 
iV = 1 state upstream of the decelerator [Sj . In the CaF 
source about 1% of all the molecules are in the state of 
interest, the |4, 0) state, and we do nothing to increase 
this fraction. 

The Stark decelerator has 100 stages and has the struc- 
ture shown in Fig.[l|^a) , the dimensions being ro = 1 mm, 
R = 1.5 mm and L = 6 mm. For the LiH experiments 
the centre of the first stage lies 135 mm from the valve 
nozzle, while for the CaF experiments this distance was 
126 mm. Four high voltage switches are used to switch 
the electrodes between two voltages, ±Vhi and ±Vio. In 
Fig.jlja) Vio — 0, but in the experiments we can choose a 
different value in order to increase the minimum electric 
field in the machine, allowing the probability of nonadi- 
abatic transitions to be controlled. 

After passing through the decelerator the time-of-flight 
profile of the molecular pulse is recorded, with 10 /is 
resolution, by cw laser induced fluorescence. The LiH 
molecules are detected 785 mm from the nozzle by ex- 
citing the A^J:+{v' = 4) - X^Y,+ {v" = 0) R(l) tran- 
sition at 367.2 nm. The CaF molecules are detected 
810 mm from the nozzle by exciting the A^ni/2(w' = 0) - 
X'^T,+ {v" 0) Q(4) transition at 606.3 nm. The fluores- 
cence is imaged onto the photocathode of a photomulti- 
plier tube operated in photon counting mode. The detec- 



tion is done at zero electric field and so cannot distinguish 
between the various M sublevels. For CaF, the ground 
state hyperfine structure is resolved and for all data we 
detect on the F" — 5 hyperfine component. Of the 11 
sublevels of this component, 3 correlate to Mjv = at 
high field, while the remaining 8 are distributed equally 
amongst the other 4 Mjy states (as shown in Fig.lTl. 



B. Results 

Figure [8] compares experimental results for the deceler- 
ation of LiH with the results of simulations. In this exper- 
iment the electrodes are switched between V^i = 9.5 kV 
and ground, and the decelerator is operated at a syn- 
chronous phase angle of 51°. The beam entering the de- 
celerator has a mean speed of 437 m/s and a temperature 
of 1.5 K. The figure shows the LiH signal as a function 
of time, the zero of time being the moment when the ab- 
lation laser fires. The curves are normalized to the peak 
of the corresponding curves obtained when the voltage is 
applied constantly to the decelerator. The lower trace in 
the figure shows the experimental data, which was also 
discussed in 6 . The narrow peak that arrives at ap- 
proximately 2.5 ms corresponds to molecules decelerated 
from 420 m/s to 213 m/s, while the broad distribution 
centred near 1.7 ms is due to the molecules that are out- 
side the phase-space acceptance of the decelerator. The 
upper trace shows the result of a numerical simulation 
of the experiment where we assume that the molecules 
always follow the electric field adiabatically. The longi- 
tudinal velocity distribution used in the simulations is 
a Gaussian distribution with centre and width set equal 
to those measured in the experiments. The simulation 
does a good job of predicting the arrival time, width and 
amplitude of the decelerated bunch of molecules. The 
shape of the non-decelerated distribution is also repro- 
duced satisfactorily, but its amplitude is not. Instead, 
the simulation predicts that this amplitude should be 
approximately 2.5 times larger than observed in the ex- 
periment. 

Next, we modify the trajectory simulations to include 
the effects of nonadiabatic transitions in a simple way. 
The transition probability is first calculated as a func- 
tion of the longitudinal position in the decelerator us- 
ing the methods described in Sec.|TT] for the voltage, rise 
time and fall time used in the experiment. We record 
the position of every molecule every time the decelera- 
tor switches, and then look up the probability of each 
molecule remaining in the |1,0) state after each switch. 
Multiplying together all these probabilities we obtain the 
total probability for each molecule to reach the end of the 
decelerator without switching state. We then take the 
detection probability to be proportional to this survival 
probability. The middle trace in Fig. |8] shows the result 
of the modified simulation. As expected from the discus- 
sion above, we find that nonadiabatic transitions lead to 
loss of the non-decelerated part of the beam, but have no 
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FIG. 8: (Color online) Time of flight profiles of LiH molecules 
decelerated from an initial speed of 420 m/s to a final speed 
of 213 m/s. The deceleration parameters are Viii = 9.5 kV, 
Vio = 0, and (j)o — 51°. (a) Lower trace, experimental data, 
(b) Middle trace, simulation including nonadiabatic loss, (c) 
Top trace, simulation for the adiabatic case. In all cases the 
signal is normalized to that obtained in dc mode. The traces 
are offset for clarity. 



effect on the decelerated bunch of molecules. The modi- 
fied simulation matches the experimental data very well 
indicating that this is the right loss mechanism, and that 
the simple procedure we have applied is a good one in 
this case. This procedure involves a number of simplifi- 
cations. We take the transition probability to be inde- 
pendent of the transverse coordinates. We have neglected 
the hyperfine structure, which is a good approximation 
for LiH as discussed above. By simply multiplying tran- 
sition probabilities, we are assuming that the superpo- 
sition state rapidly decoheres, which is reasonable given 
the field inhomogeneity. Finally, instead of calculating 
the trajectories of molecules that have changed state we 
assume that they are lost rapidly from the decelerator 
by being pulled onto the electrodes. In principle, a LiH 
molecule transferred to the |1, 1) state during one switch 
could be transferred back to the |1, 0) state during a sub- 
sequent switch and so be recaptured in the decelerator 
before being lost. However, a LiH molecule in the |1, 1) 
state is lost from the decelerator on the timescale of a 
few switch periods, and since the transition probability 
is quite small everywhere, there is very little chance of 
this recapture process occurring. 

Figure [9]^ a) shows experimental time-of-flight data ob- 
tained for CaF in the Af — 4: state, when the deceler- 
ator high voltage is Vhi ~ 14 kV. Profile (i) is the one 
obtained when the high voltage is applied statically to 
all the electrodes of the decelerator. The molecules are 
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FIG. 9: (Color online) (a) Experimental time of fiight profiles 
for CaF in the A/" = 4 state with a synchronous phase angle 
of 0° (no deceleration), a synchronous speed of 370 m/s, and 
an applied voltage of Vhi = 14 kV. (i) No switching, (ii) Vio = 
1.45 kV. (in) Vio = 0. (iv) Vio = 0.5 kV. (b) Simulations of 
these experiments. All the results share a common baseline. 



focussed through the decelerator and the time-of-flight 
profile is a direct reflection of the velocity distribution in 
the beam. For this data, the mean speed is 368m/s and 
the translational temperature is 5K. The amplitudes of 
all the profiles are normalized to the peak of this one. 

For the other three profiles in Fig.|9f^a) the decelera- 
tor is switched with a synchronous phase angle of 0° (no 
deceleration) and a synchronous speed of 370 m/s. The 
only parameter that changes between these three profiles 
is the value of V\o. For profile (ii) the decelerator was 
switched between Vhi — 14 kV and Vio — 1-45 kV. In this 
case, the minimum electric field on the axis of the deceler- 
ator is raised to 12.0kV/cm, and the corresponding min- 
imum splitting between the M = and M = ±1 states is 
2.0 X 10* rad/s. This is larger then the maximum rotation 
rate (see Fig.lsl) and we expect the evolution to be ap- 
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proximately adiabatic in this case. We observe a narrow 
central peak that corresponds to a bunch of molecules 
that are close in phase space to the synchronous molecule. 
This bunch is trapped in the moving potential well of the 
decelerator. We also see four smaller narrow peaks, two 
on either side of the main peak, and these correspond 
to bunches of molecules with speeds close to the syn- 
chronous speed but positions that are one and two po- 
tential wells ahead or behind the central bunch. Further 
out again are a pair of broader peaks corresponding to 
molecules that are not phase-stable, meaning that they 
are not trapped in the potential wells of the decelerator. 
Nevertheless, they are focussed through the decelerator, 
though not as effectively as in the case of static fields 
because the position-averaged field is only half as great 
when the decelerator is switching. We note that increas- 
ing V\o beyond this value of 1.45 kV does not change the 
time-of-flight profile, indicating that the fuUy-adiabatic 
regime has been reached. For profile (iii) the decelerator 
was switched between 14 kV and ground, i.e. Vio — 0. 
For this case, the minimum field is 2.7kV/cm and the 
minimum splitting is 1.0 x 10^ rad/s. The number of non 
phase-stable molecules is drastically reduced as we would 
expect since the maximum rotation rate is considerably 
larger than the minimum Stark splitting and so the loss 
due to nonadiabatic transitions is severe. The number 
of phase-stable molecules is also reduced but by a much 
smaller factor. This is what we expect for a synchronous 
phase angle of zero because, taking snapshots at each of 
the switch moments, the phase-stable molecules follow 
closed trajectories centred on = in the phase-space 
plot of Fig.[l]jc), and only some fraction of those trajecto- 
ries overlap the zone around —90° where transitions are 
likely. The total number of molecules transmitted by the 
decelerator in this case is 22% of the number transmitted 
in the adiabatic case. For profile (iv) the value of l^o is 
500 V. In this case, the minimum field is 5.3kV/cm and 
the minimum Stark splitting is 4.0 x 10^ rad/s. This is 
still not sufficient to avoid transitions and we see that 
the degree of loss is intermediate between cases (ii) and 
(iii) as we would expect. Again, the broad wings of the 
time-of-flight profile are reduced in amplitude far more 
than the central narrow peak. 

Figure ^h) shows the results of simulating these CaF 
experiments. We simulate the motion of molecules in all 
9 Mjv substates of = 4, and then use our knowledge of 
the correlation between the Mn states at high field and 
the F states at low field to calculate how many molecules 
are detected in _F = 5 (see Fig.[7|. Profile (i) is the result 
expected when the field is static, and serves to normalize 
the other profiles, exactly as for the experimental data. 
Almost all of the molecules transmitted are those in the 
Mtv = and 1 states, as expected since these are the 
low-field seeking states. Profile (ii) is the result obtained 
when the decelerator is switching and we assume that 
molecules always follow the field adiabatically. Again, 
although we simulate all the M^v states only the and 1 
states make a substantial contribution to the final signal. 



The result obtained agrees very well with profile (ii) of 
Fig.|9][a) , showing that the molecules do indeed follow the 
field adiabatically once Vio ~ 1.45 kV, and showing that 
their motion through the decelerator is well understood 
in this case. The central peak in the simulation is slightly 
broader and slightly smaller than in the experiments, but 
we find these features to be sensitive to the initial spatial 
distribution of molecules in the source, which is quite 
uncertain. 

To obtain profile (iii) we have simulated the case of 
Vio — 0, where the loss due to nonadiabatic transitions 
is most severe. If we do this simulation in the same way 
as we did for the LiH case, assuming that molecules are 
immediately lost if they are switched out of the Mjv = 
state, we find that the only molecules transmitted arc the 
small bunch that are phase-stable. All of the un-trapped 
molecules arc lost because the transition probability is 
high over a substantial region of the decelerator. This 
is not what happens in the experiments. Molecules are 
focussed in both the Mn = and — 1 states, and 
those that enter the decelerator in a high-field seeking 
state, or are switched into a high-field seeking state, may 
subsequently be switched to a low-field seeking state and 
thereby be transmitted. Indeed, the transition proba- 
bilities are so high that we can expect many of the un- 
trapped molecules to undergo several changes of state 
as they propagate through the machine. For the simu- 
lation that leads to profile (iii), we have tried to cap- 
ture this complexity. We first calculate the position- 
dependent probability for transitions between each and 
every M state, using the same parameters as in the exper- 
iment, neglecting the dependence of the probabilities on 
the transverse coordinates, and neglecting the hyperfine 
structure. For each molecule, we assign an initial M state 
at random, solve the motion up to the moment of the first 
switch, re-assign the M state according to the molecule's 
position and the pre-calculated probabilities, solve the 
motion up to the moment of the second switch, and con- 
tinue with this procedure until the molecule crashes or ex- 
its the end of the decelerator. Unlike the adiabatic case, 
we find that molecules in all initial M states contribute 
to the final signal, though and ±1 still contribute the 
most. The simulated profile obtained this way agrees 
remarkably well with the experimental data, suggesting 
that the calculated transition probabilities are approxi- 
mately correct despite the neglect of hyperfine structure, 
that the motion of the molecules in all relevant quantum 
states is well understood, and that the approximations 
we have made are reasonable. 



IV. CONCLUSIONS 

In this paper we have shown, both theoretically and 
experimentally, that nonadiabatic transitions occurring 
in a Stark decelerator can result in a significant loss of 
molecules. When the decelerator switches, the rotation 
of the electric field vector drives transitions to other M 
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states, some of which are high-field seeking. The transi- 
tion probabihty depends on how the rotation rate of the 
field vector compares to the splitting between neighbour- 
ing M states. In order to ensure adiabatic behaviour, 
the electric field in the decelerator should always be 
large enough that the minimum splitting between states 
greatly exceeds the maximum rotation rate. This can be 
achieved by applying a bias electric field, and choosing 
switching rise and fall times that arc as slow as possible 
while still being compatible with efficient deceleration. 
For a typical deceleration experiment, the maximum ro- 
tation rate is about 10^ rad/s. This sets the scale for the 
size of the bias field required, which will depend on the 
molecular species. A magnetic field could also be used 
to suppress nonadiabatic transitions for some molecules, 
but its magnitude would need to be at least 0.01 T for 
it to be effective. This woiild not be useful for S states 
where the only relevant magnetic interaction is with the 
electron spin (if there is one). In this case it is the com- 
bination of the Zceman interaction and the spin-rotation 
interaction that removes the degeneracy at zero electric 
field, and the induced splittings are small because the 
spin-rotation interaction is small. 

We have shown that the transition probability is a 
strong function of the position in the decelerator, and 
that the decelerated group tends to be immune to nona- 
diabatic transitions because these molecules are always 
in a large field when the decelerator switches. This is 
particularly true when the synchronous phase angle is 
large. The fact that this loss mechanism distinguishes 
between the trappc;(l and untrapped molecules could be 
used to filter out only those molecules that are being de- 
celerated. This might be useful in experiments with slow 
molecules where the non-decelerated molecules present 
an unwanted background. Our experiments on the decel- 
eration of LiH and CaF show clear evidence for loss due 
to nonadiabatic transitions. As expected, they show that 
the loss is far greater for the non-decelerated group of 
molecules and that this loss can be eliminated by apply- 
ing a bias field of an appropriate size. The experimental 
results are in excellent agreement with simulations that 
include the transitions between the M states. 

Wc have presented a simple formalism for calculating 
the M-changing transition probabilities. The formalism 
is based on linear rigid rotor molecules, but with some 
modifications could also be applied to other molecules 
of interest, e.g. asymmetric rotors. Hyperfine structure 
complicates the calculation of the transition probabili- 
ties. However, if the hyperfine structure is small com- 
pared to the Stark shift at the minimum field encoun- 
tered, as for LiH for example, it can be neglected. Our 
consideration of the hyperfine structure of CaF suggests 
that although this structure is large enough to be impor- 
tant, and though we have not done the complete calcu- 
lations, the full results will not differ greatly from those 
obtained with hyperfine structure omitted. 

For molecules that are not S states, we expect A- 
doubling to suppress the probability of nonadiabatic 



transitions. Taking OH in the ^Hz/2;J = 3/2 state 
as an example, and neglecting hyperfine structure, A- 
doubling separates the low-field seeking components from 
the high-field seeking ones by almost 10^° rad/s at zero 
field. Therefore, there are never any transitions between 
the two A-doubled states. Within the upper A-doubled 
manifold, the states differing in M have very different 
Stark shifts, one being far larger than the other, and so 
transitions between these Af states, which are degenerate 
at zero field, could still result in loss. However, because 
of the relatively small splitting between the A-doubled 
components, the Stark shift becomes linear in relatively 
small fields. As a consequence, the Stark splitting of the 
two low-field seeking states is large enough to avoid nona- 
diabatic transitions for all electric fields typically found 
in the decelerator, even when there is no bias field. For 
example, this splitting is 2.6 x 10^ rad/s at 2kV/cm. 

Another interesting and relevant example is the 
I J, K) = |1, 1) state of ammonia. Its Stark shift changes 
from a quadratic to a linear dependence on the field mag- 
nitude once the Stark shift exceeds the inversion split- 
ting. In ^^NHa the inversion splitting is large and the 
Stark shift is quadratic for fields up to ~ 20kV/cm. At 
a field of 2kV/cm, the splitting between neighbouring 
M states is 1.5 x 10^ rad/s. In "ND3 on the other hand, 
the inversion splitting is far smaller and the Stark shift 
becomes linear for fields greater than ~ 5kV/cm. At a 
field of 2kV/cm, the splitting between neighbouring M 
states is 1.9 x 10^ rad/s, more than 10 times larger than 
for NH3. These observations suggest that, unless a bias 
field is used. Stark deceleration experiments on NH3 will 
exhibit loss of the non-decelerated group due to nonadi- 
abatic transitions, whereas for ND3 no such loss should 
occur. 
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Appendix A: Time evolution in an adiabatic bsisis 

Consider some general Hamiltonian which we write as 
a sum of two parts H{t) = Hi{t) + H^it), and where the 
instantaneous eigenvalues and eigenvectors of Hi{t) are 
€m{t) and \m{t)): Hi{t)\m{t)) = em{t)\m{t)) . Using a 
basis formed by these instantaneous eigenvectors, we can 
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write the Schrodinger equation for the state vector \^p) as Thus, for m! ^ m we get 



d 



+ H2H{m\i;) - *;i^(m'|-(|m)(m|^)) 

m m 

— Ht-^^ ( ("^'I ^ I'Tt) (m|?/') + {m' \m) — {ml-tpjj . (Al) 



{m'\^^\m) = (e,„ - e,n'){m'\^\m). (A6) 



Using the notation c„i = ("T'IV') obtain the following 
equation for the time evolution of the coefficients Cm- 
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^(m'|7?i + H2 - ih—\m)cm = ih 
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(A2) 



Using this result in Eq. ( A2 ) we obtain 



Taking the derivative of Hi — X^m" ^m"\'m"){m"\ with 
respect to time we get 
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Finally, writing c, 
we arrive at 



ih dHi \ ,dCm' 
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It follows that for m! ^ m 



d{m'\ 
dt 



I to). (A4) 



The last term in this equation can be put into a more 
convenient form by noting that 

|((to'|to» = = (^) \m) + (m'|||TO). (A5) 



i dHi 



ih- 



ddm' 



(A8) 

This way of expressing the time evolution of the state 
vector makes it clear what operator is responsible for in- 
ducing transitions between the instantaneous eigenstates, 
which of the eigenstates are coupled by the time- varying 
Hamiltonian, and what condition must hold for the evo- 
lution to be adiabatic. 
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